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The amount of heat transferred to green coffee beans is essential in the coffee roasting. During this pro¬ 
cess, several parameters can be used as indicators to determine the degree of roasting (color, aroma, vol¬ 
ume, bean temperature). Consequently, two predictive models, using artificial neural networks, are 
proposed to determine the quality of coffee roasting. A bean brightness model took account of bean tem¬ 
perature (simulated by a physical model) and roasting time. A second model predicting the bean surface 
area, focused on roasting air temperature and roasting time. The color changes affecting coffee beans dur¬ 
ing the process were studied experimentally in a pilot roaster equipped with a CCD video camera and a 
lighting system consisting of two small optic fiber spotlights. Arabica green coffee beans of Colombian 
origin were roasted using different air temperatures (190, 200, 210,... 300 °C), for 10 min. Two separate 
feedforward networks with one hidden layer were used to brightness and surface kinetics. The best fit¬ 
ting training data set was obtained using three neurons in the hidden layer, which enable prediction of 
brightness and bean surface kinetics with an accuracy that was at least as good as the experimental error, 
over the entire experimental range. Using the validation data set, simulations and experimental data 
were in good agreement ( R 2 > 0.98). This study showed that real-time simulations were possible so that 
the roasting process could be stopped when the simulated brightness was similar to a target fixed by the 
roast master. The model thus developed could contribute to the on-line estimation of product quality, 
thus providing a parameter to control coffee roasting. 

© 2008 Elsevier Ltd. All rights reserved. 


1. Introduction 

Coffee roasting is an unitary operation of considerable impor¬ 
tance to developing the specific organoleptic properties (flavor, ar¬ 
oma and color) which underlie the quality of coffee and guarantee 
a satisfactory beverage. However, this process is highly complex: 
during coffee roasting, moisture loss and chemical reactions (oxi¬ 
dation, reduction, hydrolysis, polymerization, decarboxylation 
and numerous other chemical changes), as well as major altera¬ 
tions (color, volume (swell), weight, bean pop, pH, volatile compo¬ 
nents and C0 2 production, etc.) occur. Finally, after these 
considerable changes, the beans must be cooled rapidly to halt 
these reactions (using water or air as a cooling agent) in order to 
prevent excessive roast that would affect product quality (Sch- 
wartzberg, 2000; Illy and Viani, 1998; Nagaraju et al., 1997; Raemy 
and Lambelet, 1982; Raemy, 1981; Dutra et al., 2001; Singh et al., 
1997; Sivetz and Desrosier, 1979; Guyot, 1993). 

In the food industry, the quality of roasted coffee is evaluated 
off-line using different parameters (mainly aroma, flavor, color, 
pH, bean pop, mass loss and volume). However, in such an indus¬ 
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trial setting, it is very difficult to estimate these parameters on¬ 
line, and in most cases the roast master has an essential role to 
play. The operating conditions are determined by the roast master 
based on off-line measurements obtained on the previous batch 
(color, aroma, flavor, etc.), and then adjusts the roasting parame¬ 
ters (air temperature and duration of the process) for the next 
batch accordingly (Hernandez-Perez, 2002). 

On the other hand, in a context of handmade roasting, the oper¬ 
ator monitors several parameters of the process; the most impor¬ 
tant being bean brightness, the moment at which the beans 
begin to crunch and finally bean volume. Roasting times are not 
standardized as in industry, but rather adjusted on-line. Today’s 
food industry suffers from a lack of accurate on-line sensors, which 
would constitute a major breakthrough in term of enabling indus¬ 
try automation. Direct on-line brightness measurements are very 
difficult in an industrial roaster, but it is possible to determine 
one or more variables which are correlated with this quality, and 
then to build a model that can predict the quality of the coffee as 
a function of process parameters. 

It is possible to built a semi-empirical model (heat and mass 
transfer model from physical laws considering chemical reactions) 
to predict bean temperature and moisture content in coffee roast¬ 
ing as a function of process parameters (Schwartzberg, 2000; 
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Hernandez et al., 2007; Heyd et al., 2007) and then use such semi- 
empirical model to predict brightness as a function of the time- 
temperature history. However, the chemical reactions leading to 
bean browning at high temperatures are extreme complex and at 
present poorly understood. Another way to model quality as a 
function of process variables is to employ an empirical model (Arti¬ 
ficial Neural Network). When modeling such a process, the ability 
of artificial neural networks (ANN) to integrate complex relation¬ 
ships between process parameters and product quality is of great 
interest. The major advantage of such networks is that they can 
model without any assumptions as to the nature of underlying 
mechanisms, and they can take account of non-linearities and 
interactions between different variables (Bishop, 1994). An out¬ 
standing feature of neural networks is their capacity to learn the 
solution to a problem from a set of examples, and provide a smooth 
and reasonable interpolation for new data. In the field of food pro¬ 
cess engineering, neural networks have recently been applied to 
correlating color to moisture content in cooked beef (Qiao et al., 
2007; Chaoxin et al., 2006). 

The aim of this work was to simulate the brightness and surface 
area of coffee beans during roasting using artificial neural net¬ 
works, in order to build a smart sensor for on-line quality measure¬ 
ment. The determination of surface area is also important when 
calculating the heat transfer coefficient between beans and air. 

2. Materials and methods 

2.1. Coffee roasting 

Colombia green coffee beans ( Arabica ) were roasted using a hot 
air flow as the heating medium. According to experimental infor¬ 
mation reported by Schwartzberg (2000); Hernandez-Perez 
(2002), experiments were carried out at a constant air velocity of 
4 y, at different fixed air temperatures (190, 200, 210, 220, 230, 
240, 250, 260, 270, 280, 290 and 300 °C) for 10 min. All experi¬ 
ments were performed in triplicate. 

A static roaster (SERVATHIN Series SV02 7817) was used to per¬ 
form roasting experiments. A diagram of the roaster is shown in 
Fig. 2. In order to simulate industrial roasters, where convection 
is the major mode of heat transfer, the beans were placed on a 
mesh grid to maintain them in a static suspension. Type K thermo¬ 
couples allowed to measure bean internal temperature and to con¬ 
trol the air temperature with a precision of ±0.5 °C. Beans were 
drilled in order to insert a thermocouple 4 mm below the internal 
temperature (see Fig. 1). Air temperature was measured at 20 mm 
of the bean temperature. Beans weight was automatically mea¬ 
sured using an electronic scale METTLER (SB16001) and controlled 
by the computer (a program written in bash (Bourne, 2002) making 



the weight acquisition each minute as it is reported by Hernandez 
et al. (2007)). During coffee roasting, air temperature, bean tem¬ 
perature and bean weight were measured (Hernandez-Perez, 
2002; Hernandez et al., 2007), and color pictures of the beans were 
also acquired on-line (Hernandez et al., 2008). 

2.2. Data acquisition system 

Fig. 2 shows the experimental system which was developed to 
monitor the brightness and surface area (projected area) of coffee 
beans during roasting. This experimental system was made up of 
standard devices used for image analysis: 

• An illumination system: a light source with two small optic fiber 
spotlights (Illuminate Technology Inc. USA), 

• An image sensor: a CCD (Charge-Coupled Device) video camera 
(SONY XC-711P) operated with a 50 mm </>25.5 objective, 

• A digitization system: a computer, equipped with a Hauppauge 
WinTV video acquisition card and a bttv 878 converter for image 
digitization. These commercial systems were deliberately used 
because of their robustness and lower cost, and particularly 
the availability of source codes for the software pilots for the 
bttv 878 component, enabling modifications for specific use, 
including process controls. 

Air temperature data acquisition was ensured as follows: 

• Thermocouples (type K) to measure air temperature with an 
accuracy of ±0.5 °C, 

• An Arcom SCB7 thermocouple conditioner connected to the per¬ 
sonal computer and, 

• An Arcom PCAD12/16H A/D converter for the acquisition of air 
temperatures in the roaster. 

Data acquisition and I/O port programs were written in C (Law¬ 
yer, 2000; Photis, 1999). On-line image data processing was per¬ 
formed using a program written in Octave (Eaton, 2001) and an 
algorithm managed under bash (Bourne, 2002). 

As shown in Fig. 2, the video camera was installed outside the 
roasting cell and events were visualized through a borosilicate 
glass window (the image distortion induced was negligible). Bttv- 
grab software was used for image acquisition (Walter, 2001 ). This 
software has the advantage that it can be used with a program 
written in bash , enabling simple on-line image processing and pro¬ 
cess control. Images were stored on the hard disk in a ppm (porta¬ 
ble pixel map) format at regular intervals (one image every twenty 
seconds) with a definition of 480 • 640 pixels. The image analysis 
system comprised three stages: acquisition (visualization of 
objects, numeration), processing and then extraction of the 
information. It is important to note that all measurements (air 
temperature, brightness and surface) were acquired on-line, so 
that process decisions could be taken real-time. 

2.2.1. Image processing 

The aim of image processing was to enable on-line control the 
roasting process based on quality parameters (brightness and sur¬ 
face), and to determine the degree of roasting. In order to simulate 
industrial conditions as closely as possible, it was chosen not to 
control the light level. The images acquired thus depended on 
the experimental conditions, resulting in heterogeneous images 
(each point of the image not receiving either the same quantity 
or quality of light). For this reason, before each experiment, the 
system was calibrated using a color palette which was photo¬ 
graphed under exactly the same conditions as the beans (Hernan¬ 
dez-Perez, 2002). The brightness of each color was measured using 
a Minolta CR-200 colorimeter. According to Vrhel and Trussel 
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(1999) and Schettini et al. (1995), a neural network can be used to 
determine the look up table (LUT) between Red-Green-Blue coor¬ 
dinates given by the acquisition system, and the normalized 
brightness given by the colorimeter (Hernandez-Perez, 2002). 
Thresholding made it possible to isolate the beans from the back¬ 
ground, so that the mean brightness of beans could be calculated. 


3. Experimental results 

3.1. Brightness 

Fig. 3 shows the brightness (hr) values measured using the CCD 
video camera and neural network LUT during different coffee 
roasting experiments. The behavior was similar at all tempera¬ 
tures, there being a rapid reduction between 20 and 60 s, followed 
by an almost symmetrical growth (around 60-100 s) and then an 
uninterrupted exponential decline. 

It was possible to break this behavior down into four different 
stages: 

(1) bean color remained the same during an initial period (per¬ 
iod from 0 to 20 s), 

(2) when the internal temperature of beans reached 100 °C (see 
Fig. 4), their color darkened slightly (after approximately 
20-60 s), which could be attributed to the vaporization of 
unbound water, 

(3) above 160°C, the beans were markedly lighter (period 
between 60 and 100 s), 

(4) the Maillard and pyrolytic reaction then began, resulting in 
gradually darker beans that finally reached the color of 
burned coffee (determined visually). 


Under high air temperatures (^260°C) brightness changes 
were more rapid during the process. Brightness curves (with a 
R = 0.9999) thus confirmed that air temperature and time are 
two crucial factors during the roasting process. 

3.2. Projected surface areas 

Distances on pictures were counted in numbers of pixels. These 
depended on the experimental conditions and particularly on the 
focal length of the objective and the distance between objects 
and the camera. In order to calibrate surface area measurements, 
a picture was acquired which contains a black disk (30 mm) on a 
white background placed instead of the coffee beans (600 mm 
from the camera objective). The measured dimensions of the disk 
were not significantly affected by its position on the picture, be¬ 
cause of the long focal length of the objective. Calibration showed 
that a pixel corresponded to 0.0125 mm 1 2 3 4 of the disk. The image 
sensor had a fixed aperture and focus for all experiments. 

After calibration of the pixel surface, the bean projected surface 
area was measured from the number of pixels belonging to the 
beans. Fig. 5 shows the bean surface area as a function of air tem¬ 
perature and roasting time. The increase in surface area was be¬ 
tween 15% and 70% at air temperature of between 190 and 
300 °C. Schwartzberg (2000) found a volume expansion of be¬ 
tween 50% and 120% (corresponding to respectively 30% and 70% 
of the projected surface area increase if swelling was considered 
as isotropic) at air temperatures of between 270 and 550 °C, using 
a different roasting system. At higher air temperatures (over 
270 °C), the increase in surface area remained stable, but at that 
point, optimum bean brightness was exceeded. Coste (1968) spec¬ 
ified that the volume of coffee beans no longer increased above a 
temperature of 280 °C. Our experimental findings showed (Fig. 5) 
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Fig. 3. Brightness experimental kinetics for three different air temperatures and their repetitions. 



that coffee beans only began to swell when their internal temper¬ 
ature reached 100 °C (see Fig. 4 for comparison), suggesting that 
swelling is mainly a consequence of internal steam pressure in 
the beans. 

4. Modeling 

Some semi-empirical models based on the laws of chemical 
kinetics have been proposed to model brightness as a function of 
process parameters (Broyart et al., 1998; Krokida et al., 2001). 
However, these models were not appropriate to fit our data, and 
especially brightness kinetics which exhibit a complicated non- 
monotonous behavior (Fig. 3). Rather than constructing a semi- 
empirical model combining different physical laws that involve 
chemical reactions, it is turned to using a mathematical approach 
involving neural networks. However, these empirical models have 
a narrower validity range, but require only a limited number of 



Fig. 5. Experimental kinetics of bean surface area at different air temperatures 
(190-300 °C). 


simple arithmetic operations for simulation, and can be easily 
incorporated in control software (Hernandez-Perez et al., 2004). 

According to previous results (Hernandez-Perez, 2002), the use 
of different architectures as input variables in the input layer had 
been proposed to predict the brightness and surface area kinetics, 
the input layer was considered from different parameters combi¬ 
nations (bean and air temperatures, time, moisture content, rate 
of exothermic heat production) (data not shown). The best results 
for the brightness were obtained based on two input variables: 
roasting time (t) and bean temperature T b , simulated using the 
dynamic model proposed by Hernandez et al. (2007), because 
Schwartzberg (2000) also suggested that bean temperature is a 
good indicator of the level of coffee roasting (see Fig. 6a). It was 
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also possible to use the experimental bean temperature and roast¬ 
ing time as input variables to model brightness, but this solution is 
not realistic in an industrial setting because it is very difficult to 
obtain this parameter in an industrial roaster. In order to surface 
area kinetics, the best results were computed simply from the air 
temperature T a and roasting time t as input variables (see 
Fig. 6b). These surface curves were easy for computed because 
the curves exhibit a simple behavior when compared with bright¬ 
ness. The surface variable is an important parameter to compute 
bean temperature, so that it was not possible to use bean temper¬ 
ature as input variable. 

The artificial neural networks were inspired by the study of 
neurosciences that today has applications in the food industry 
(Qiao et al., 2007; Chaoxin et al., 2006; Broyart et al., 1998) and 
notably in the specific area of image analysis for the real-time def¬ 
inition of product quality (Boillereaux et al., 2007; Park and Chen, 
2000). Neural networks are capable of learning the dynamics of a 
process from experimental data, taking account of non-linearities 
of the system and correlations between variables. As natural neu¬ 
rons, they are determined to a great extent by connections be¬ 
tween two elements, with each connection between two neurons 
having a coefficient (weight). This notion of weighting enables 
modulation of the sign transmitted between two neurons as a 
function of the state of the link computer synaptic that links them 
(Hagan and Menhaj, 1994). 

Neurons are organized in several layers that are interconnected 
to a given architecture. These models use classical neural networks 
of multi-layer perceptron formed by three elements, this being typ¬ 
ical for the approximation of functions. These elements are formed 
by the input layer, hidden layer and output layer. Each element of 
the layer is connected to each neuron input through the weight 
matrix. The optimum architecture is determined from the minimal 
value of root mean squart error obtained (Demuth and Beale, 
1998). 

In order to calculate the stimulation Sj of a neuron in the hidden 
layer, it is necessary to consider the activation A, of each neuron of 
the input layer, which is then multiplied by its corresponding 
weight Py. The bias Bj is then added to regulate the threshold of 
activation of the neuron (Demuth and Beale, 1998). 

Sj = jr(Pij-A i ) + Bj (1) 

i= 1 

The activation function is then applied to calculate Ap 



Experimental database 

Fig. 7. Experimental versus simulated brightness values for all the learning 
database. 


A = -- 2 - „ N - 1 (2) 

3 1 + exp(-2 • Sj) y 3 

These calculations can be expressed more simply under the matrix 
form (Dornier et al., 1998) 


S c -P c A e + B c (3) 

4c =fc(Sc) (4) 

S s =P s A c + B s (5) 

4 S =fs(S s ) (6) 


where A e , is the input standardized by the model; S c , is the stim¬ 
ulation of the hidden layer,/ c , is the activation function (hyperbolic 
tangent sigmoid transfer-function) (Eq. (2)), A c , is the activation of 
the hidden layer and S s , is the stimulation of the output layer,/ s , is 



Fig. 6. Architecture of the neural networks proposed to determine the product quality using roasting time as an input variable. 
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Table 1 

Statistical parameters obtained (weights and biases) in the best model for the 
brightness RMSE leaming = 0.7392, RMSE testing = 0.739 


Weights 

Pi 


-1.6133 


-7.4940 



6.7830 


3.4317 



7.5821 


4.5031 

Po 

Biases 

0.2230 


-1.9935 

1.9021 

B ! 




5.3858 

-3.0258 

-3.7821 

b 2 




-0.1192 


the activation function (linear transfer-function), to A s =f s . The 
most commonly used transfer functions for backpropagation are 
the logarithmic sigmoid (logsig), hyperbolic tangent sigmoid 
(tansig) and linear (purelin) (Demuth and Beale, 1998). But other 
different transfer function can be created and used with backprop¬ 
agation if desired. According to the literature (Demuth and Beale, 
1998; Hernandez-Perez et al., 2004), the authors proposed that 
the appropriate notation is used in the two-layer (hyperbolic tan¬ 
gent sigmoid transfer-function/ linear transfer-function) network 
shown in the Eqs. (2) and (6). In the last-layer, the linear trans¬ 
fer-function is used because the outputs of the network are unlim¬ 
ited (network outputs can take on any value). In addition, in 
previous results (Hernandez-Perez, 2002), the notation tansig/ 
purelin transfer function was better than the logsig/purelin trans¬ 
fer function. 

The coefficients of the network (weight P and bias B), the num¬ 
ber of neurons in the hidden layer and the number of iterations of 
the algorithm of optimization are calculated during the training 
stage, thus minimizing a root mean square error (RMSE) between 
simulated and experimental data in the learning base. The opti¬ 
mum model is that which introduces minimal error. During the 
present study, the neural networks toolbox of the Matlab software 
was used (Demuth and Beale, 1998) with a Levenberg-Marquardt 
optimization algorithm, considered by Hagan and Menhaj (1994) 
as being the most efficient, since this algorithm was designed to 
approach second-order training speed without having to compute 



Experimental brightness 



Fig. 9. Experimental versus simulated surface area kinetics for all the learning 
database. 


the Hessian matrix. To test the robustness and then the prediction 
ability of the models, the experimental database was split into the 
learning and test databases. Experimental data (brightness and 
surface kinetics) were measured at 12 different air temperatures 
(190-300 °C) in triplicate, thus generating 36 kinetics for bright¬ 
ness and 36 kinetics for bean surface area. These experimental 
database were split into learning (67%) and testing (33%). 24 kinet¬ 
ics (air temperatures of 190, 210, 220, 240, 250, 270, 280 and 
300 °C) comprised the learning database, enabling the calculation 
of optimal weights and biases, while the remaining 12 kinetics 



Time (s) 


Fig. 8. Testing of the neural model regarding brightness, (a) experimental and simulated brightness values for all the test database, (b) brightness values as a function of time 
at air temperature of 200, 230, 260 and 290 °C: (*) experimental data and (-) simulated data. 
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Table 2 

Statistical parameters obtained (weights and biases) in the best model for the bean 
surface RMSE\ eaming = 0.1598, RMSE testing = 0.159 


Weights (P f ,P 0 ) 

Pi 

4.5578 


-0.2101 


-3.4513 


-4.1150 


5.4978 


2.3551 

Po 

-1.9052 

0.3703 

0.8104 

Biases (B\,B 2 ) 

B\ 



-5.0779 

5.5003 

-5.5690 

b 2 



-1.4239 


(at 200, 230, 260 and 290 °C), comprised the test database, en¬ 
abling validation of the model and testing its capacities for general 
implementation. 

4A. Brightness model 

A neural network with 3 neurons in the hidden layer (involving 
13 coefficients: 9 weights and 4 bias) was found to be efficient in 
predicting the brightness of coffee beans during the roasting. 
Fig. 7 presents a comparison between experimental and simulated 
brightness values using the learning database. The neural network 
was well fitted to the non-monotonous behavior of the data during 
all phase of the process. The correlogram shown in Fig. 7 shows 
that the residuals were small for all experiments in the learning 
base (R 2 > 0.98). 

Table 1 shows the RMSE for learning and testing respectively, as 
well as the obtained statistical parameters (Pi,P 0 ,B^and B 2 ) of the 
proposed model. Consequently, these coefficients are used in the 
artificial neural network model to simulate brightness values. 

The model thus calculated was then validated using the kinetics 
in the test database (fresh data). Fig. 8a compares experimental 
and simulated brightness values for roasted beans in the test data¬ 
base (roasting air temperatures of 200, 230, 260 and 290 °C). The 
neural model also fitted these unknown data well, including the 
initial phases of the curves (first 100 seconds) which are complex 


a 0 10 20 30 40 50 60 70 



and very difficult to predict using a model including assumptions 
regarding physical-chemical mechanisms (Fig. 8b). The residuals 
in the test database (Fig. 8a) were small and their distribution 
well-balanced. The coefficient of regression ( R 2 = 0.98) confirmed 
this satisfactory agreement. The model thus developed was able 
to predict the brightness of beans as a function of process param¬ 
eters throughout the entire time-temperature experimental do¬ 
main. The similarity of the coefficients of regression for the 
learning and test database confirmed the predictive capacity of this 
model. 

4.2. Model for bean projected surface area 

In order to predict bean surface area kinetics, the best network 
was similar to that described above (three neurons in the hidden 
layer). Fig. 9 shows the experimental and simulated surface area 
values for the learning database. The data fitted well ( R 2 = 0.99 
for the learning database). 

Table 2 shows the RMSE for learning and testing respectively, 
as well as the obtained statistical parameters (Pi,P 0 ,Bi and B 2 ) 
of the proposed model. Consequently, these coefficients are used 
in the artificial neural network model to simulate bean surface 
values. 

In the same way as for brightness. Fig. 10a compares experi¬ 
mental and simulated surface area values for the test database 
(fresh data). The model exhibited a satisfactory ability for predic¬ 
tion (R 2 = 0.99 for the test database). Fig. 10b presents the evolu¬ 
tion of experimental and simulated surface area for kinetics 
included in the test database (200, 230, 260 and 290 °C). The ability 
to simulate an unknown situation was confirmed by the compara¬ 
ble coefficients of determination. 

5. Conclusion 

This study proposes two neural network models that can pre¬ 
dict brightness and bean surface areas during coffee roasting as a 
function of process parameters, and exhibit a good ability for 
generalization. 


b 0 100 200 300 400 500 600 



Time (s) 


Fig. 10. Testing of the neural model to reveal increases in bean surface area during roasting, (a) experimental data as a function of the values simulated for all the test 
database, (b) surface area as a function of time at air temperature of 200, 230, 260 and 290 °C: (*) experimental data (with repetitions) and (-) simulated data. 
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The model to predict bean surface area may be of assistance 
better understanding heat and mass transfer during drying, be¬ 
cause the contact surface between beans and heat media, which 
is a key parameter when modeling product temperature, can be 
computed simply from the projected surface area (insofar as swell¬ 
ing is roughly isotropic). 

Brightness is a crucial factor when determining coffee quality 
during the roasting process, and the model proposed enables the 
on-line computation of this parameter, based on the roasting time 
and bean temperature as calculated using a dynamic semi-empir¬ 
ical model that takes account of physical laws (Heyd et al., 2007; 
Hernandez et al., 2007). 

The models thus developed may enable the implementation of 
smart sensors for on-line quality determination in industrial 
plants: 

• the exchange surface between bean and air will be computed 
from the roasting and air temperature using the proposed neural 
network. 

• bean temperature can be computed from the roasting time, air 
temperature and exchange surface area using the dynamic heat 
and mass transfer model, 

• the validity of the bean temperature thus computed will be con¬ 
firmed by comparing the measured and computed outlet roast¬ 
ing air temperature (a simple measurement in an industrial 
setting), 

• the computation of product brightness using the neural network 
is then possible and make it possible to halt the process when 
the target quality has been attained. 

It will also be interesting to take in account evolutions in bright¬ 
ness during cooling so that the quality of the final product can be 
optimized. This will be made possible by computing the end of 
the roasting phase in order to achieve the desired quality at the 
end of cooling. 

References 

Bishop, C.M., 1994. Neural networks and their applications. Review Science 
Instrument 65, 1803-1832. 

Boillereaux, L., Cadet, C., Le Bail, A., 2007. Thermal properties estimation during 
thawing via real-time neural network learning. Journal of Food Engineering (1), 
17-23. 

Bourne, S., 2002. Bash (free software). In Bourne Again SHell. Free Software 
Foundation [ On-line ]. <http: //www.gnu.org/software/bash/bash.html>, 

consulted June 2008. 

Broyart, B., Trystram, G., Duquenoy, A., 1998. Predicting colour kinetics during 
cracker baking. Journal of Food Engineering 35, 351-368. 

Chaoxin, Z., Da-Wen, S., Liyun, Z., 2006. Correlating colour to moisture content of 
large cooked beef joints by computer vision. Journal of Food Engineering (4), 
858-863. 

Coste, R., 1968. Le cafeier. G.P. Maisonneuve et Larose, XIV ed., 11, rue Victor-Cousin, 
Paris (V e ). 


Demuth, H., Beale, M., 1998. Neural Network Toolbox for Matlab, User’s Guide 
Version 3. The MathWorks Inc., USA. 

Dornier, M., Heyd, B., Danzart, M., 1998. Evaluation of the simplex method for 
training simple multilayer neural networks. Neural Computing and 
Applications 7, 107-114. 

Dutra, E.R., Oliveira, L.S., Franca, A.S., Ferrez, V.P., Afonso, R.J., 2001. A preliminary 
study on the feasibility of using the composition of coffee roasting exhaust gas 
for the determination of the degree of roast. Journal of Food Engineering 47, 
241-246. 

Eaton, J.W., 2001. Octave (free software). In GNU Octave, University of Wisconsin 
[On-line], <http://www.gnu.org/software/octave/, consulted june 2008. 

Guyot, G., 1993. Torrefaction: mecanismes, transformation physiques et chimiques. 
In Journees du CAFE, CIRAD-CP a Montpellier. France. 

Hagan, M., Menhaj, M., 1994. Training feedforward networks with the Marquardt 
algorithm. IEEE Transactions on Neural Networks 5 (6), 989-993. 

Hernandez, J.A., Heyd, B., Irles, C., Valdovinos, B., Trystram, G., 2007. Analysis of the 
heat and mass transfer during coffee batch roasting. Journal of Food Engineering 
78,1141-1148. 

Hernandez, J.A., Heyd, B., Trystram, G., 2008. On-line assessment of brightness and 
surface kinetics during coffee roasting. Journal of Food Engineering 87, 314- 
322 . 

Hernandez-Perez, J.A., 2002. Etude de la torrefaction: modelisation et 

determination du degre de torrefaction du cafe en temps reel. Ph.D. Thesis in 
Ecole Nationale Superieure des Industries Agricoles et Alimentaires. France. 

Heyd, B., Broyart, B., Hernandez, J.A., Valdovino-Tijerino, B., Trystram, G., 2007. 
Physical model for the heat and mass transfer in the coffee roasting. Drying 
Technology 61 (6-7), 212-222. 

Hernandez-Perez, J.A., Garda-Alvarado, M.A., Trystram, G., Heyd, B., 2004. Neural 
networks for the heat and mass transfer prediction during drying of cassava and 
mango. Innovative Food Science and Emerging Technologies 5 (1), 57-74. 

Illy, A., Viani, R., 1998. Espresso Coffee. Academic Press, San Diego, CA 92101, US. 

Krokida, M.K., Oreopoulou, Z.B., Marinos, D., 2001. Colour changes during deep fat 
frying. Journal of Food Engineering 48, 219-225. 

Lawyer, S.D., 2000. Serial-HOWTO. In the Linux Serial HOWTO [On-line]. <http:// 
www.linuxdoc.org/HOWTO/Serial-HOWTO.html>, consulted January 2005. 

Nagaraju, V.D., Murthy, C.T., Ramalakshmi, K., Srinivasa, R.P.N., 1997. Studies on 
roasting of coffee beans in a spouted bed. Journal of Food Engineering 31, 263- 
270. 

Park, B., Chen, Y., 2000. Real-time dual-wavelength image processing for poultry 
safety inspection. Journal of Food Process Engineering 23, 329-351. 

Photis, M.G., 1999. Coffee-HOWTO.txt. In software coffee-HOWTO, University 
Cambrige [On-line]. <http://www.faqs.org/docs/Linux-mini/Coffee.html># 
ssl.l, consulted January 2007. 

Qiao, J., Ngadi, M.o., Wang, N., Gariepy, C., Prasher, S.O., 2007. Pork quality and 
marbling level assessment using a hyperspectral imaging system. Journal of 
Food Engineering (1), 10-16. 

Raemy, A., 1981. Differential thermal analysis and heat flow calorimetry of coffee 
and chicory products. Thermochimica Acta 43, 229-236. 

Raemy, A., Lambelet, P., 1982. A calorimetric study of self-heating in coffee and 
chicory. Journal of Food Technology 17, 451-460. 

Schettini, R., Barolo, B., Boldrin, E., 1995. Colorimetric calibration of color scanners 
by back-propagation. Pattern Recognition Letters 16, 1051-1056. 

Singh, P., Singh, R., Bhamidipati, S., Singh, S., Barone, P., 1997. Thermophysical 
properties of fresh and roasted coffee powders. Journal of Food Process 
Engineering 20, 31-50. 

Sivetz, M., Desrosier, N.W., 1979. In: Coffee Technology. AVI Publishing Co., 
Westport, pp. 226-264. 

Schwartzberg, H.G., 2000. Modelling bean heating during batch roasting of coffee 
beans. In: Welti-Chanes, J., Barbosa-Canovas, G., Aguilera, J.M. (Eds.), 
Engineering and Food for the 21st Century. CRC Press, LLC, London, New York, 
Boca Raton. 

Vrhel, M.J., Trussel, H.J., 1999. Color scanner calibration via neural network. In: IEEE 
ICASSP, vol. 6, pp. 3465-3468. 

Walter, J., 2001. bttvgrab-0.15.10.tar.gz (free software). In bttvgrab [On-line]. 
<http://www.garni.ch/bttvgrab/>, consulted January 2007. 


